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Abstract 

A general Lorentz-invariant parameterization for the long-range part of the 0uf3f3 
decay rate is derived. Combined with the short range part this general parameter- 
ization in terms of effective B — L violating couplings will allow it to extract the 
Qv(3f3 limits on arbitrary lepton number violating theories. Several new nuclear ma- 
trix elements appear in the general formalism compared to the standard neutrino 
mass mechanism. Some of these new matrix elements have never been considered 
before and are calculated within pn-QRPA. Using these, limits on lepton number 
violating parameters are derived from experimental data on 76 Ge. 
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1 Introduction 



Double beta decay has been proven to be one of the most powerful tools to con- 
strain B — L violating physics beyond the standard model [1] . In recent years, 
besides the most restrictive limit on the effective neutrino Majorana mass [1,2], 
stringent constraints on several theories beyond the Standard Model such as 
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.R-parity violating [3-5] as well as conserving [6] SUSY, leptoquarks [8], left- 
right symmetric models [9] and compositeness [10,11] have been derived (for 
a review see [1]). 

While the neutrino mass limit is based on the well-known mechanism exchang- 
ing a massive Majorana neutrino between two standard model V — A vertices, 
the effective vertices appearing in the new contributions involve non-standard 
currents such as scalar, pseudoscalar and tensor currents. 

Thus we felt motivated to consider the neutrinoless double beta decay rate in 
a general framework, parameterizing the new physics contributions in terms of 
all effective low-energy currents allowed by Lorentz-invariance. Such an ansatz 
allows one to separate the nuclear physics part of double beta decay from 
the underlying particle physics model, and derive limits on arbitrary lepton 
number violating theories. The first step of this work, treating the long-range 
part, is presented here. Although the general decay rate is independent of 
the underlying nuclear physics model, to extract quantitative limits values 
for nuclear matrix elements are needed. First limits are derived using matrix 
elements calculated in proton-neutron (pn) QRPA, partially already available 
in the literature and partially calculated here for the first time. 



2 General Formalism 



We consider the long-range part of neutrinoless double beta decay with two 
vertices, which are pointlike at the Fermi scale, and exchange of a light neu- 
trino in between. The general Lagrangian can be written in terms of effective 
couplings ejg, which correspond to the pointlike vertices at the Fermi scale so 
that Fierz rearrangement is applicable: 

£ = %^4-^ + £^4} (i) 

with the combinations of hadronic and leptonic Lorentz currents .P a = uO a d 
respectively jp = eOpv of defined helicity. The operators O a ^ are defined as 

0y-A = 7 M (l-7s) 
O v+ A = 7 M (l+7 5 ) 

Os-p = (1 - 7s) 

O s+ p = (1+75) 

<^ = ^,7,](l-75) 

Or„ = ^[7M.7j(l + 75). (2) 
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The prime indicates the sum runs over all contractions allowed by Lorentz- 
invariance, except for a = (3 = V — A. Note that all currents have been scaled 
relative to the strength of the ordinary (V — A) interaction. 



The effective Lagrangian given in eq. (1) represents the most general low- 
energy 4-fermion charged- current interaction allowed by Lorentz invariance. 
The interpretation of the effective couplings e^, however, depend on the spe- 
cific particle physics model. Nevertheless one realizes the following general 
feature. Using only the SM fermion fields^] and working in the Majorana ba- 
sis for the neutrinos (u :— vl + z/£) it is easily seen that all currents involving 
operators proportional to (1 + 75) violate lepton number by two units, i.e. 
the corresponding must also be lepton-number violating. (Such LNV are 
easily found, an example is given by R-parity violating supersymmetry treated 
in ref. [3-5].) 

The double beta decay amplitude is proportional to the time-ordered product 
of two effective Lagrangians (see Fig. 1): 



Q2 

T(£ ( i)£ {2 )) = -YT{j v _ A J ] v _ A j v „ A J ] v _ A 

+€^p4jv-a4-a + 44ip J ihJ\}- (3) 

While the first term (contribution (a) in Fig. 1) corresponds to the standard 
model (SM) like neutrino exchange, and the 3rd term (contribution (c) in Fig. 
1), which is quadratic in e can be neglected, only the 2nd term (contribution 
(b)) is phenomenologically interesting. For this term one has to consider two 
general cases: 

1 ) The leptonic SM V — A current meets a left-handed non SM current jp 
with (3 = S — P,T L . For this contribution the neutrino propagator is 

~ ft^^ft = -P^s (4) 

with the usual left- and right-handed projectors Pl/r = -^p-- This expres- 
sion is proportional to the unknown neutrino Majorana mass m u <~ 0.5eV, 
for which no lower bound exists. Therefore no limits on the corresponding 
parameters can be derived. 



' The use of only SM neutrinos does not really imply a loss of generality here, see 
the discussion on V + A currents below. 
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2) The leptonic SM V — A current meets a right-handed non SM current jp 
with (5 — S + P, V + A, T R . For this contribution the neutrino propagator is 



+ m 





(5) 



which is proportional to the neutrino momentum ~ p F ~ lOOMeV with 
the nuclear Fermi momentum p F , and thus will produce stringent limits on 
corresponding e^. 

Taking these considerations into account, we are left with three interesting 
contributions discussed in the following section. With the present half-life limit 
of the Heidelberg-Moscow experiment T^ 13 > 1.2 • 10 25 y [1] and considering 
only one at a time (evaluation "on axis") 



where G k denotes the phase space factors given in [12] and \ME\ the nuclear 
matrix elements discussed in the following. Note that evaluating "on axis", 
compared to the arbitrary evaluation, neglects interference terms of the differ- 
ent contributions. Although this is expected to be a small effect (see below), 
this remains to be discussed in the next step. 



3 Calculational details and limits 

3. 1 SM meets j v+A J v+A and 3v+aJ ] V -a 

This contribution has been considered already in the context of left-right 
symmetric models [12,15,9]. For sake of completeness we repeat the updated 
results of [15] here in our notation: ty%A < 7 - 9 • 10 " 7 > ^v-A < 4 - 9 • 10 ~ 9 for 
the full calculation with arbitrary evaluation. Using s-wave approximation 
and "on axis" evaluation as was done in this work the limits are reduced to 
^vXa < 7.0 • 10~ 7 , ey^A < 4-4 • 10~ 9 . This example confirms the expectation 
that these assumptions will only slightly (less than 10 %) affect the result. 

It is worthwhile discussing the following subtlety in interpretation when com- 
paring our ansatz and the work of [12]. Doi et al. [12] calculate the decay 
amplitude, writing down the Lagrangian of an explicitly left-right symmet- 
ric model. Therefore, their calculations do not only treat SM neutrinos, but 
contain also right-handed neutrinos. Nevertheless, the derivation of the decay 
amplitude for Qvf3(3 decay is the same in both calculations, only some care 
is required when going from one notation to the other. For example, in the 



{ei) 2 G 0k \ME 



(6) 
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lepton- number violating (A), defined in [12] as 

(\) = X&ejVej (7) 
j 

where A ~ (mw L /mw R ) 2 , LNV is due to the product of mixing matrices 
J2'jU e jV e j, which vanishes identically if LNV goes to zero. Thus, our (LNV) 
e v+A corresponds to (A) of Doi et al. [12], (and not to A). This example shows, 
that the use of only SM neutrinos in the derivation does not imply a loss of 
generality in our results, if the source of LNV of the particle physics model 
under consideration is carefully identified. 



3.2 SM meets js+pJg+p and js+pJg-p 



Using s-wave approximation for the outgoing electrons and some assumptions 
according to [13,4] one gets 



,(3) 
P 

Rm e G A 

The phase space factor is defined 



ME%Xl = -ME§+J = ~^§-{M r + \m gt ). (8) 



5(ei + e 2 + M f - M i )de 1 de 2 (9) 



and the matrix elements (summation over nucleons a, b is suppressed) are 

M GT> = (o;iM^KVlo l + > (10) 

M T > = (0f\h T/ {(a a f ab )((j b f ab ) 

-\{^ b )}r:r^t). (11) 

hp and h T > are neutrino potentials defined as 

, 2 i? 2 r 4 j (gr ab )/V) 

h R = / dqq p=r— , 12 

7rmpi uj(uj + E) 



2 i? 2 ? 2 / 2 (g 2 ) 2 



nmpJ uj(uj + E) 
Q 



-3^ 3l (qr ab )}. (13) 

' ab 
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Here R denotes the nuclear radius, mp the proton mass, q and pi are elec- 
tron energies and momenta, F (Z,ei) and is the Fermi function. Further uo = 

\Jq 2 + m 2 ,, q = \q\, r = f/r and jk(qr) are spherical Bessel functions, (ou + 
E)^ 1 is the energy denominator of the perturbation theory. The form factors 
F?(0) = F?(q 2 )/f(q 2 ) with f(q 2 ) = (1 + q 2 /m 2 A )- 2 (m 2 A = 0.85 GeV) have 
been calculated in the MIT bag model in [14], G A — 1-26 and G v — 1- 

Inserting the numerical value of the matrix elements M.qt' an d At T > (see [4] 
and Tab. 1), one derives ef+p, 4-p < 11 ' 10 ~ 8 - 



3.3 SM meets jT R Jr R an d 3t r Jt l 

In the tensor part the decay rate depends on the phase space Goi an d new 
matrix elements not considered in the literature. 

For the hadronic T R contribution one gets under the assumptions used above 
ME^l = -a x -M GT > + a x M T > (14) 



with 

AT{ 3 \0)Gv(l-2m P (G w /G v )) 
ai ~ G\R^ e • (15) 

Again T-p^O) has been taken from [14] and the strength of the induced weak 
magnetism (Gw/Gy) = — is obtained by the CVC hypothesis. The 

involved nuclear matrix elements have been calculated in the QRPA-approach 
of [15,16]. Inserting the values obtained for the special case of 76 Ge (see Tab. 
1) yields eg < 1.7 • 1CT 9 . 

For the hadronic Tl contribution in leading order of (1/m.p) one finds 

ME^l = a 2 M F > - az(M T » + \m gt ») (16) 



with 



M F > = <0;|/^t+t+|0+> (17) 

M GT » = (0+|^((7>l)r a + r+|0+) (18) 

M T n = (tf\uh T »{{<7 a r ah )((7 h r ah ) 

4(*a* 6 )}r^+|0t). (19) 
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The neutrino potentials are 



nmpJ^ (co + E) 



Hrp" 



2 R 3 f , 2 /V) f 2 ■ , , 



m P J (u + E) 
-3— ii(?r ab )}. (21) 



and 



4(2f 2 (3) (0)-T 1 (3) (0))G v 

" 2 = cp^ ' (22) 

47f ) (0)(Gp/^) 

tt3 = G^R%^ • (23) 

Here nas been taken from [14], Gp/GU = 2mp/m^ and uo is the neu- 

trino energy. Again the matrix elements have been calculated in the model of 
[15,16]. A limit of e^P < 7.3 • 10~ 10 has been obtained. Factors R have been 
arbitrarily absorbed into the definition M. G t" an d M. T " to get dimensionless 
quantities. One should notice that the corresponding factor included in a 3 
compensates this choice. 



4 Conclusion 



We have presented a general parameterization for the long range part of the 
neutrinoless double beta decay rate in terms of effective couplings. The result- 
ing bounds are summarized in Tab. 2. Combined with the short range part 
and contributions of derivative couplings, this parameterization will give the 
double beta decay constraints for arbitrary lepton number violating theories 
beyond the SM. The next step should include these contributions and discuss 
interference terms. 
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(a) (b) 




Fig. 1. Feynman graphs of the general double beta rate: The contributions a) - c) 
correspond to the long range part, the contribution d) is the short range part (to be 
discussed elsewhere). 



M GT' 


2.95 


M F > 


-0.663 


M GT " 


8.78 


M T > 


0.224 


M T " 


1.33 


Table 1 



Relevant Nuclear Matrix Elements calculated in pn-QRPA 



e V+A 
e V-A 


4.4 • 10~ 9 


V+A 
e V+A 


7.0 • 10~ 7 


£ S+P 
t S-P 


1.1 • 10~ 8 


e S+P 
e S+P 


1.1 • 10" 8 


JTR 
e TL 


6.4- 10~ 10 


JTR 
e TR 


1.7- 10~ 9 



Table 2 

Limits on effective B — L violating couplings evaluated "on axis 



9 



